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Abstract. In order to examine the excitation and ioniza- 
tion mechanism of extraplanar diffuse ionized gas (DIG) 
we have obtained optical longslit spectra of seven edge- 
on spiral galaxies. In four objects the brightest emission 
lines can be traced out to distances of typically 1.5 kpc 
above the disk. For NGC 1963 and NGC3044 hne ra- 
tios such as [Nil] A6583/Ha or [Sii] A6717/Hq: as well as 
[O III] A5007/H/3 could be measured for the halo DIG. This 
allows us to discuss the DIG in the halo of these objects 
in the framework of diagnostic diagrams. For these two 
objects, the line ratios of [O ill] A5007/II/3 decrease with 
increasing 1 2 ] , different from the recently reported trend in 
NGC 891 (Rand |l998l ). We find that emission lines from 
the DIG in the disks are in good agreement with pho- 
toionization models using a dilute radiation field. How- 
ever, with increasing \z\ these models fail to predict the 
measured [O i] A6300/Ha and He i A5876/Hq! line ratios for 
NGC 1963. Diagnostic diagrams reveal for NGC 1963 the 
need for a second ionization mechanism of the halo DIG 
(besides photoionization) . This additional source could be 
shock ionization. The same diagrams demonstrate an in- 
termediate classification for NGC 3044. Plots of [S 11] /[N 11] 
vs. emission measure reveal significant changes towards 
the halo and seem to trace local small scale density fluc- 
tuations of the extraplanar DIG. 
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1. Introduction 

Over the last decade a widespread and extended diffuse 
component of ionized hydrogen contributing 25 - 50 % to 
the total Ha flux has been found well outside of classical 
Hii regions in several face-on and edge-on galaxies (e.g., 
Lehnert & Heckman 1995, Rand 1996| , Greenawalt et al. 



1997, or Dettmar 199S). Although there is significant evi- 



dence that the presence of the extraplanar diffuse ionized 



gas (DIG, sometimes also called warm ionized medium, 
WIM) is correlate d with ongoing star formation in the 
disk (e.g., Dettmar 1992), important aspects of its ioniza- 
tion and excitation are not well understood. While of all 
conventional ionizing sources only UV radiation of mas- 
sive stars provides sufficient power to keep these thick 
gas layers ionized (including the Reynolds-layer of the 
Milky Way (Reynolds 1995 )) the resulting photoioniza- 
tion models have to explain the long mean free paths for 
the UV photons and the observed emission line ratios. 
However, the relatively strong emission observed in the 
[N 11] AA 6548, 6584 lines (Dettmar & Schulz |1992| , GoUa et 
al. 1996, Rand 1997) cause a problem for pure photoioniza- 
tion models and an additional heating contribution, e.g. 
from the kinetic energy provided by supernova, is likely 
(Dahlem et al. 1997). An even more puzzling observation 
was recently presented for the DIG in NGC 891 which is 
up to now the best studied galaxy in this regard (Rand 



1998). In complete contradiction to photoionization mod- 



els here [O ill] /H/3 is rising with distance above the galac- 
tic plane for \z\ > 1 kpc. Diagnostic diagrams for emis- 
sion lines from DIG in outfiows of dwarf galaxies (Mar- 
tin 1997) also demonstrate the need for ionization and/or 



excitation processes beyond pure photoionization. In this 
paper we report results from spectroscopic observations 
of a small sample of edge-on galaxies. These new spectra 
cover a larger wavelengths range extending earlier work 
to the blue range of the optical spectrum. This allows us 
to discuss the results in the framework of diagnostic di- 
agrams. The paper is structured as follows: In Sect. 2-3 
some essential information concerning observations and 
data reduction strategies, the determination of line ra- 
tios, and the use of diagnostic diagrams is given. Sect. 4 
shows representative results which are discussed briefly for 
each galaxy with respect to the photoionization models by 
Mathi s ( |1986| , hereafter Ma86) and Domgorgen & Mathis 
(1994, hereafter DM94). Finally, our findings are sum- 
marized in Sect. 5 with special emphasis to the trend of 
[S 11] / [N 11] , and the proposed dependence of line ratio vari- 
ations with changing halo metallicities. 
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2. Observations and data reduction 

The dataset consists of spectra of the normal late- type spi- 
rals NGC 1963, IC 2531, NGC 3044, NGC 4302, NGC 4402, 
NGC4634, and of the irregular galaxy NGC 2188. AU ob- 
jects in the sample are seen nearly edge-on (see Table |l|). 
Slit positions were determined from narrow band Ha im- 
ages on which diffuse Hli gas at hig h galactic latitudes 
was detected (Rossa & Dettmar 2000 ). 

The spectroscopic data were obtained during Jan- 
uary/February 1995 with the ESO 1.52m telescope in com- 
bination with the Boiler & Chivens spectrograph at La 
Silla (Chile). A Ford Aerospace CCD (2048 x 2048 pix- 
el) with a pixel size of 15 x 15 and a spatial resolution 
of 0.82" pix~^ was used during the observations. Grating 
23 provided a dispersion of 126Amm~^ which resulted 
in a coverage of a wavelength range from 3500 A -7450 A. 
We determined the spectral resolution to 4.6 A by mea- 
suring the FWHM of an unbroadened emission line in a 
HeAr-calibration spectrum. The total integration time per 
slit position was 90 minutes. 

Data reduction was performed using the IRAF and 
MIDAS software packages. All data were bias subtracted 
and dome flats were used to reduce small scale sensitiv- 
ity variations. An extinction correction was performed by 
adopting the atmospheric extinction curve for La Silla 
(Tiig 1977 ) . Spectroscopy of the standard star Hiltner 600 
yielded the wavelength dependent response function. We 
divided all object frames by this function to correct for 
wavelength dependent sensitivity variations. Subsequently 
a flux calibration of our spectra was applied. 

The slit length of 4f2 allowed a precise night sky sub- 
traction. This reduction process did not influence the 
emission lines of HeiA5876 and [Oi]A6300 since they are 
sufficiently redshifted from the atmospheric lines for our 
objects. Due to the large wavelength range we could not 
apply a single continuum fit to the whole spectrum. There- 
fore the continuum level was estimated and subtracted row 
by row from the total emission of each galaxy using lin- 
ear fits to the continuum on each side of prominent emis- 
sion lines. The data for N GC 21 88 are already presented 
by Domgorgen & Dettmar ( 1997 ) and will not be further 
discussed here. As a consistency check these spectra were 
reanalysed and the previous results were fully reproduced. 

3. Line ratios and diagnostic diagrams 

For all galaxies listed in Table 1 emission lines are traced 
down to S/N = 2a. Detected spectral features for all galax- 
ies are [N ii] AA6548, 6583, Ha and [S ii] AA6717, 6732. Line 
ratios are measured using a block average algorithm. The 
block size was fixed by the achieved S/N (at least 2) to 
reach simultaneously the best possible spatial resolution. 
For the disk and the halo of a galaxy we have chosen block 
sizes of 1-2 and 3-4 rows, respectively. 

For the galaxies NGC 1963 and NGC 3044 additional 
emission lines could be detected from the halo DIG, 



namely [O ii] A37271], H/3, [O ill] AA4959, 5007, HeiA5876, 
and [O i] A6300. Thus the extraplanar DIG in these ob- 
jects can be discussed for the first time in the framework 
of diagnostic diagrams. A typical spectrum of the DIG in 
NGC 1963 is given in Fig. 2. 

A comparison with photoionization models of Ma86 
and DM94 can clarify some properties of the sources 
which keep the DIG ionized. Only the line ratios of 
[On] A3727/[Oiii] A5007 and [Oiii] A5007/H/3 have to be 
corrected for interstellar extinction. The precise values are 
obtained following Osterbrock ( 1989| ). Unreddened line ra- 
tios can be computed by assuming an average interstellar 
extinction curve and adopting H/? as nebular reference 
emission line. The extinction can be determined from the 
measured Balmer decrement. It is assumed that such a 
correction is negligible for the remaining line ratios be- 
cause only pairs of emission lines in the red spectral region 
had been chosen which are located in a narrow wavelength 
region (< 200 A) so that differential reddening is not sig- 
nificant. 

It is noteworthy that an extinction correction does not 
influence a clear separation of different excitation mech- 
anisms significantly (Baldwin et al. p_981 ). The dered- 
dening for [O ii] A3727/[0 III] A5007 and [O ill] A5007/H/3 
amounts to 11% for the disk emission and 4% for the halo 
emission. For the remaining line ratios the reddening cor- 
rection is less than 1% (disk and halo). 

In extragalactic objects basically four different ion- 
ization mechanisms can be clearly distinguished. Firstly 
"normal" photoionization by O stars which is present 
mainly in Hii regions (star forming regions). Secondly 
shock ionization present in objects which are excited by 
supernova remnants (SNR) or by superwinds. Finally we 
can also distinguish between photoionization by extremely 
hot stars such as the central stars of planetary nebulae 
and objects which are photoionized by a power-law con- 
tinuum source (AGN). To probe the ionization mechan- 
sism(s) of the extraplanar DIG it is useful to plot line ra- 



tios by pairs on a logarithmic scale (Baldwin et al. 1981 
Veilleux & Osterbrock 1987). Since the above mentioned 



mechanisms produce characteristically different spectra, 
all objects cover different locations when plotted into a 
single diagnostic diagram. 

In order to verify our findings obtained from line ratios 
and diagnostic diagrams we examined also the broaden- 
ing of line widths at extraplanar latitudes. Therefore we 
plotted a kinematical diagnostic diagram (line ratio vs. 
line width (FWHM)) which indicates the presence of ad- 
ditional excitation mechanisms, such as shocks or stellar 
winds. 

4. Results and discussion 



^ Due to low spectral resolution the [O ll]-doublet cannot 
be separated. We therefore assume an average wavelength of 
3727 A. 
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Table 1. Sample overview 



Galaxy 


R.A. 


Dec. 


Type 


D 


z" 


inclination 


rriR 


DIG morphology® 




(J 2000) 


(J 2000) 




[Mpc]" 




i^ 


[mag]'^ 




NGC 1963 


05'' 33" 12.8" 


-36°23'59" 


Sc 


17.7 


0.004410 


85° 


12.11 


diffuse 


NGC 3044 


09''53'"39.8" 


+Q1°34'46" 


SBb 


17.2 


0.004310 


84° 


11.65 


bright, diffuse 


IC2531 


09''59'"55.7" 


-29°36'55" 


Sb 


33.0 


0.008460 


90° 


11.54 


1 filament 


NGC 4302 


12''21'"42.4° 


+14°36'Q5" 


Sc 


18.8" 


0.003690 


88°" 


11.83" 


faint, diffuse 


NGC 4402 


12''26'"07.9" 


+13°06'46" 


Sb 


22.0" 


0.000790 


74°" 


12.09" 


diffuse 


NGC 4634 


12''42'"40.4" 


+14°17'47" 


Sc 


19.1" 


0.000390 


83°" 


11.91" 


bright, diffuse, extraplanar 


















+ filaments 



for Ho = 75kms"^Mpc"\ ^ BottineUi et al. 1984, " Teerikorpi et al 



1992 



The Surface Photometry Catalogue of the ESO-Uppsala Galaxies; Lauberts A. and Valentijn E.A. 
" NASA Extragalactic Database 

^ Third Reference Catalogue of Bright Galaxies; de Vaucouleurs et al., 
" Rossa& Dettmar 2000 



1989 



1991 




5000 
wavelength [A] 



Fig. 2. NGC1963: Representative intensity plot (z = 3", 
corresponding to 260 pc) for slit position si. Rest wave- 
lengths (in [A]) of relevant emission lines are marked. This 
spectrum shows the line and continuum emission. 



4.1. NGC1963 

NGC 1963 is a Sc galaxy with an inclination angle of 85° 
(edge-on) at a distance of 17.7 Mpc (see Table ||). The 
corresponding slit positions si and s2 are overplotted on 
the Ha image and are displayed in Fig. |l|. Some diffuse 
emission near the nucleus in the disk-halo interface and 
bright extended Hii regions in the disk are visible. The 
origin of the z-scale was determined from the continuum 
emission of the disk. 

Although generally plotted as a function of z, we also 
present line ratios as a function of emission measure (EM) . 
This provides additional information about the depen- 



dence of line ratios and the mean squared densities of free 
electrons in the emitting region. Line ratios for positions si 
and s2 are plotted in Figs. |~|. Table | gives an overview 
of averaged line ratios for the disk, the northern, and the 
southern halo. Representative logarithmic diagnostic dia- 
grams for slit position si are displayed in Fig. ^. Diffuse 
Ha emission can be traced spectroscopically in the north- 
ern halo out to 1.4 kpc and out to 1.3 kpc in the south. The 
line ratios [N ii] /Ha and [S ii] /Ha measured near the mid- 
plane of NGC 1963 are typical for Hii regions and show 
no conspicuous features. Both line ratios increase towards 
the outer halo (|z| > 500 pc). They are well fitted by the 
DM94 model taking into account a composite model with 
20% (X.dge = 0.95) + 80% (Xedge = 0.1). X^dge is the frac- 
tion of neutral hydrogen at the outer border of the model 
nebula. If X^dge is high (radiation bounded case) nearly 
all emitted energy is absorbed by the geometry and only 
few stellar photons ionize hydrogen at the edge. A lower 
value of A"edge (matter bounded case) indicates that a sig- 
nificant fraction of ionizing photons escapes the geometry. 

At both slit positions the averaged value for Hei/Ha 
is 0.05 ± 0.01 (cf. Table ||). With respect to the error bars 
the DM94 model fails in predicting the correct ratios no 
matter which parameter set is chosen. This is also true 
for the shock sensitive line ratio [Oi]/Ha at position s2. 
It follows that there is most likely a second excitation 
mechanism besides photoionization. As the diagnostic di- 
agrams with [Oi]A6300 emission (see Fig. ^) reveal, this 
extraplanar source could be shock ionization. 

The line ratios of [O ill] /H/3 decrease with increasing 
z I (Fig. ^ . This pattern is easily explained in photoiniza- 
tion models (e.g., Ma86) by the ionization stratification 
expected from the ionization potentials involved. How- 
ever, recent observations in NGC 891 (Rand 1998) show 
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Fig. 3. NGC 1963: Line ratios of [S ii] A6717 / Ha, [N ll] A6583 / Ha (left) and [O i] A6300 / Ha, [O i] A6300 / [O ill] A5007, 
and Hei A5876 /Ha (right) along slit position si. 5" correspond to 430 pc. The lower left panel displays the variation 
of line ratios as a function of emission measure EMhq- Representative mean errors of [Sii]/[Nii] (halo and disk) are 
overplotted. 
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Fig. 4. NGC 1963: Same ratios of [Sii] A6717/Ha, [N ii] A6583 / Ha (left) and [O i] A6300 / Ha, 
[O i] A6300 / [O III] A5007, and HeiA5876/Ha (right) but this time for slit position s2. Plots in the lower left 
are the same as for si. Representative mean errors (halo/disk) are plotted for [Sll]/[Nll]. Note the tight correlation 
between [O i]/Ha and Hei/Ha in the lower left panels of each figure. 
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15 20 



Fig. 5. NGC 1963: Measured line ratios of [O ii] A3727 / [O ill] A5007 and [O ill] A5007 / H/3 along slit position si (left) 
and s2 (right) . 





log [01I]/[0III] 




log [01I]/[0I11] 




Fig. 6. NGC 1963: Relevant diagn ostic d iagrams for slit position si. The dotted line and filled triangles represent the 
shock model from ShuU & McKee (1979) with shock velocities between 80 -100 km s~^. Filled circles denote the disk 
data and open circles address the halo component. Crosses represent observed shock ionized objects used by Baldwin 
et al. (1981). Areas enclosed by solid lines are fitted by Ma86 and DM92. 
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Table 2. Averaged line ratios for the disk and halo of 
NGC 1963. Halo values are given for distances perpendic- 
ular to the plane as indicated below. 



slit 1 


disk 


halo (north) 


halo (south) 






z = 1.3 kpc 


z = 0.7 kpc 


[IN IIJ/ rid 


U.iO ± U.Ui 


n -J- not:: 


r\ Qd -i- r\ no 

u.oD =b u.uy 


[D IIJ/rlQ; 


n 1 o 4- n m 
u.iy ± U.UI 


n 1 n 1 o 
U.Do ± u.iy 


n oo -1- f\ no 
U.zy zt U.Uy 


rie 1 / xlQ: 


U.UO ± U.Ui 






[O i]/|OniJ 


U.i 1 ± U.Uo 






[\J l\/ rid 


U.Uo ± U.UI 






[Oni]/H/3 


2.20 ± 0.18 


1.53 ± 0.29 


0.83 ± 0.52 


[Oii]/[Oiii] 


2.06 ± 0.34 


3.89 ± 1.57 


6.57 ± 2.00 


slit 2 


disk 


halo (north) 


halo (south) 






z = 0.9 kpc 


2 — 1.0 kpc 


[Nll]/Ha 


0.18 ± 0.01 


0.26 ± 0.12 


0.53 ± 0.19 


[Sii]/Hq 


0.16 ± 0.01 


0.44 ± 0.15 


0.42 ± 0.26 


Hei/Ha 


0.05 ± 0.01 






[Oil/[Oni] 


0.13 ± 0.04 


0.33 ± 0.16 


0.15 ± 0.08 


[0 i]/Hq 


0.04 ± 0.01 


0.12 ± 0.05 


0.09 ± 0.05 


[0 iii]/H/3 


1.79 ± 0.23 


0.78 ± 0.44 


0.65 ± 0.48 


[On]/[Oni] 


2.90 ± 0.64 


6.88 ± 2.51 


6.91 ± 2.63 



an increase of the line ratio of [Oiii]/H/3 with increasing 
distance from the midplanc and therefore from the most 
likely location of the ionizing sources. 



4.2. I C 2531 

This galaxy of type Sb is perfectly edge-on and shows no 
prominent DIG emission at either of the slit positions. 
Only a single filament can be detected 20" NE of slit s2 



(cf. Rossa& Dettmar 200C) 



Table 3. Averaged line ratios for the disk and the halo 
area of IC2531. Halo values and their corresponding z- 
distance are given below. Note that the halo data for slit 1 
are only upper limits. 



slit 1 


disk 


halo (north) 


halo (south) 






z = 1.1 kpc 


z = 1.0 kpc 


[Nii]/Hq 


0.35 ± 0.17 


2.02 ± 2.02 


0.57 ± 1.54 


slit 2 


disk 


halo (north) 


halo (south) 






z = 1.0 kpc 


2 = 1.1 kpc 


[Nll]/Ha 


0.38 ± 0.07 


0.80 ± 0.33 


0.21 ± 0.16 



Because the slits do not cut any bright H ii region the 
lack of emission lines is not surprising. For both slit posi- 
tions only the hne ratios of [N ii] A6585 / Ha (Fig. 7) could 
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Fig. 7. IC2531: Line ratios of [Nil] A6583/Ha for si and 
s2 (upper panel). Additionally, plots of line ratio vs. EMho 
together with mean error bars for [Nil] /Ha are presented 
(lower panel). 5" correspond to 800 pc. 

be determined (see also Table ^) and therefore no diagnos- 
tic diagram could be obtained. 

In view of a single line ratio and relatively large er- 
ror bars (low S/N ratio) it is not reasonable to make any 
statements on possible excitation mechanisms of the DIG 
in IC2531. The obtained values for the disk region seem 
to be reproducible by DM94. It is not possible to distin- 
guish between a matter or radiation bounded geometry 
because both parameter ranges are able to fit the data. 
Ha line emission can be traced at si out to 1.4 kpc in the 
northern and out to 1.0 kpc in the southern halo. At s2 the 
detection of diffuse Ha emission reaches out to 1.3 kpc in 
the mean. 

4.3. NGC 3044 

Besides NGC 1963 this galaxy shows the strongest dif- 
fuse extraplanar emission. Again, an Ha image from 
Rossa & Dettmar was used to classify the DIG morphol- 

For NGC 3044 
and diagnostic 



ogy as "bright" and "diffuse" (cf. Fig, 
line ratios are presented in Figs. p|~|10 
diagrams are shown in Fig. |ri|. Averaged line ratios are 
given in Table ^ In this galaxy diffuse extraplanar Ha 
emission is detectable at slit si up to ± 15" which cor- 
responds to 1.3 kpc. At s2 the extraplanar DIG can be 
traced up to 1.6 kpc in the northern and up to 1.1 kpc 
in southern halo. All line ratios except that of [Sii]/Ha 
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Fig. 8. NGC 3044: Line ratios of [S ii] A6717 / Ha, [N ii] A6583 / Ha (left) and [O i] A6300 / Ha, [O i] A6300 / [O ill] A5007, 
and HeiA5876/Ha (right) along slit position si. 5" correspond to 415 pc. The lower left panel contains plots for 
different line ratios vs. EMua, including representative mean errors for [Sll] / [Nil]. 
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Fig. 9. NGC 3044: Same line ratios of [Sii] A6717/Ha, [N ii] A6583 / Ha (left) and [O i] A6300 / Ha, 
[O i] A6300 / [O III] A5007, and HelA5876/Ha (right) but this time for slit position s2. 5" correspond again to 
415 pc. Line ratios vs. EMho are also plotted for s2. Again, note the constant behaviour of Hei/Ha and [Oi]/Ha in 
the lower left panels. 
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Fig. 10. NGC 3044: Measured line ratios of [O ii] A3727 / [O ill] A5007 and [O ill] A5007 / H/3 along slit position si (left) 
and s2 (right) . 
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Fig. 11. NGC 3044: Relevant diagn ostic diagrams for slit position si. The dotted line and filled triangles are from 
shock models (ShuU & McKee 197E ) assuming shock velocities between 80 -100 km s~^. Filled circles denote the disk 
data a nd op en circles address the halo component. Crosses represent observed shock ionized objects used by Baldwin 
et al. ( 1981 ). Predicted values of Ma86 and DM92 are enclosed by solid lines. 
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Table 4. Averaged line ratios for the disk and the halo 
of NGC3044. Ratios for the halo with an asteriks (*) are 
measured within ±0.8 kpc. 



slit 1 


disk 


halo (north) 


halo (south) 






z = 1.1 kpc 


z — 1.2 kpc 


[IN IIJ/ rid 


U.zo rt U.Ui 


U.Oi ± U.iD 


U.OI ± U.Uo 


[o Ilj/rlQ; 


n 1 o -1- n f\A 
U.iy ± U.U4 


n QO 1 n 1 /I 
U.oz zt U.i4 


U.oo it U.U4 


rle I/xIq: 


U.U4 zh U.Ui 


n no ^ n oq 
U.Uo ± U.Uo 


n n^j _i_ n no 
U.Uo ± U.Uz 


[U l]/[OlllJ 


U.i ( ± U.Uo 


U.2< ± U.U9 


U.o4 ± U.io 


\r\ t1 /H^,* 
[\J IJ/rlQ: 


n no 1 n m 
U.Uz dz U.UI 


U.U/ ± U.U4 


U.Uo =t U.U4 


[Oiii]/H/3 


0.91 ± 0.10 


0.74 ± 0.29 


0.49 ± 0.38 


[Oii]/[Oiiil 


2.92 ± 0.27 


6.94 ± 2.02 


6.11 ± 1.50 


slit 2 


disk 


halo (north) 


halo (south) 






z = 1.2 kpc 


2 = 0.9 kpc 


[Nll]/Ha 


0.24 ± 0.01 


0.28 ± 0.07 


0.35 ± 0.04 


[Sii]/Hq 


0.19 ± 0.01 


0.41 ± 0.09 


0.64 ± 0.11 


Hei/Ha* 


0.04 ± 0.01 


0.07 ± 0.04 


0.08 ± 0.06 


[Oil/[Oni]* 


0.06 ± 0.03 


0.17 ± 0.04 




[0 i]/Hq* 


0.04 ± 0.01 


0.05 ± 0.03 




[0 iii]/H/3 


1.66 ± 0.10 


0.87 ± 0.27 


1.15 ± 0.37 


[On]/[Oni] 


2.45 ± 0.42 


6.39 ± 1.38 


5.73 ± 1.56 



at slit s2 for NGC 3044 are well reproduced by the pho- 
toionization model from Ma86 or DM94 assuming a mat- 
ter bounded case with Xcdgo = 0.10 and q= 0.001 {q is 
proportional to the ratio of ionizing photon density to 
electron density cubed and also proportional to the ion- 
ization parameter U). In these models ionizing sources 
are 05 stars with temperatures of T, = 4.5 ■ 10'' K. This 
means that most of the hard Lyman continuum radiation 
from the star forming regions can escape and ionize the 
medium at high galactic latitudes. The above mentioned 
line ratio [S ii]/Ha reaches values of 0.80 ± 0.16 at z = 11". 
This finding is in good agreement with the data obtained 
by Lehnert & Heckman ( 19951 ) . In diagnostic diagrams in- 
cluding the [Oi] emission line (Fig. 11) the data for the 
halo fall into the gap between H ii regions and shock ion- 
ized objects. Additionally, the separation of disk and halo 
components is not as clear as in NGC 1963. The line ratio 
of [SiiJ/Ha for the southern halo at s2, the position of 
our data in diagnostic diagrams (Fig. 11), and a moder- 
ate gradient in line ratios reveal the hybrid character of 
NGC 3044 regarding the excitation mechanism. Although 
photoionization seems to be the main ionizing source in 
this galaxy contributions due to shocks cannot be ruled 
out. Out to our detection limit we find again a decreasing 
[Oiii]/H/3 ratio with increasing \z\. 



4.4. NGC43O2 

No diffuse extraplanar line radiation is detectable in 
NGC 4302 at slit position si above 2cr. Slit s2 cuts the 
outer south western part of an extended H 11 region. There- 



fore only faint Ha and [N 11] A6583 emission from the disk 
can be detected. No empirical diagnostic diagram could be 
obtained. Averaged values for [N 11] A6583/IIa are shown 
for -4" < |z| < +4" in Table ||, with 4" corresponding 
to 360 pc. The DM94 model fails in predicting the mea- 
sured data while the "simplified" Ma86 model fits with 
respect to the error bars values up to 0.62 well. This as- 
sumes a diluted radiation field (logq = —4) and 05 stars 
(T, = 4.5 • 10** K) as ionizing sources. For values larger 
than 0.62, as observed in the outer disk, the Mathis mod- 
els fail, too. 



Table 5. Averaged line ratios for the center, the north- 
ern (N), and southern (S) outer disk of NGC 4302 at sht 
position s2. 



slit 2 


disk center 


outer disk (N) 


outer disk (S) 


[Nii]/Hq 


0.44 ± 0.06 


0.78 ± 0.24 


0.99 ± 0.24 



4.5. NGC 4402 

An R-band image of NGC 4402 reveals that this galaxy is 
rich in dust and shows some spectacular dust filaments 
emerging from the south eastern and western parts of 
the disk. The corresponding Ha image gives evidence that 
these filaments are connected with H 11 regions (star form- 
ing regions) inside the disk. Our slit position has been 
chosen such that the star in the north of the galaxy falls 
onto the slit (Fig. |l|). Despite of a "diffuse" DIG classi- 
fication no emission lines in the blue wavelength domain 
could be detected. The measured line ratios are presented 
in Fig. 12 and averaged values are given in Table O. 



Table 6. Averaged line ratios for the disk and halo area 
of NGC 4402. 



slit 1 


disk 


halo (north) 


halo (south) 






z = 1.9 kpc 


2 = 1.3 kpc 


[Nii]/Hq 


0.30 ± 0.05 


0.76 ± 0.34 


1.07 ± 0.29 


[Sn]/Ha 


0.15 ± 0.04 


0.50 ± 0.24 


0.62 ± 0.18 



For the disk region of NGC 4402 the DM94 model re- 
produces the observed ratios well. In order to fit the data 
one has to choose a matter bounded model (Xcdge = 0.10) 
and a photon field with logq = —3. Roughly 20% of the 
ionizing photons escape from the H 11 regions and are able 
to ionize sulfur and nitrogen in the halo. Higher values 
(< 0.9) for the halo are reproduced (with respect to the 
error bars) by the model of Ma86 using 05 stars as ioniz- 
ing sources and a still softer radiation field (logq — —4). 
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Fig. 12. NGC4402: Line ratios of [S ii] A6717/ Ha and 
[N ii] A6583 / Ha including variations of [Sii]/[Nii] vs. 
EMho for the only slit position si. 5" correspond to 
530 pc. Positive z- values denote the northern, and neg- 
ative the southern halo. 



Values of 0.62 for [N ii]/Ha and 0.61 for [S ii]/Ha are pre- 
dicted. Ratios exceeding 0.9 lead to a model failure. 

The steep gradient in [Nil] /Ha or [Sii]/Ha reaching 
values of 1 or higher at high \z\ cannot be explained in 
the framework of dilute photoionization models and may 
indicate again the need for an additional heating source. 

16. NGC4634 

On the Ha + [N ii] image for NGC 4634 a bright and extra- 
planar DIG layer is visible. Additionally, several filaments 
emerge mainly from the north-eastern part of the disk (see 
Fig. |l|). Here again the slits had been positioned such that 
they contain stars for accurate positioning. Slit si cuts 
the brigthest H ii region « 5" north west of an extended 
dust cloud and s2 covers an area of fainter Ha intensity 
«14" north west of si. As in NGC 4402, only the spec- 
tral Hues of [N II] AA6549, 6583, Ha, and [S ii] AA6717, 6732 
are measurable at both slit positions above 2a. The cor- 
responding line ratios are plotted in Figs. 13 and 14. Due 
to the non-detection of emission lines in the blue wave- 
length region no diagnostic diagrams could be created. In 
NGC 4634 DIG emission can be traced up to 1.2 kpc at slit 
position si and up to 0.9 kpc at s2. Measured line ratios 
of [Nil] /Ha for the disk area can be reproduced using the 
DM94 model and setting Xodgc again to 0.10 and logq to 



Fig. 13. NGC 4634: Line ratios of [S ii] A6717 / Ha and 
[Nil] A6583/Ha for slit si. 5" correspond to 465 pc. The 
lower panel displays [Sii]/[Nii] vs. EMhq in combination 
with corresponding mean errors. 

Table 7. Averaged line ratios for the disk and halo area 
of NGC 4634. 



slit 1 


disk 


halo (north) 
z = 1.2 kpc 


halo (south) 
2 = 1.2 kpc 


[Nii]/Hq 
[Sn]/Ha 


0.39 ± 0.02 
0.22 ± 0.02 


0.79 ± 0.21 
0.69 ± 0.19 


0.80 ± 0.18 
0.62 ± 0.18 


slit 2 


disk 


halo (north) 
z = 0.8 kpc 


halo (south) 
2 = 0.8 kpc 


[Nil] /Ha 
[Sn]/Ha 


0.34 ± 0.01 
0.18 ± 0.02 


0.50 ± 0.07 
0.51 ± 0.08 


0.85 ± 0.16 
0.37 ± 0.07 



—4. Only 10 % of hydrogen is neutral at the border H°, H+ 
that is 31 % of the ionizing photons can escape the geom- 
etry. The model fails again for values larger than 0.37. 
[Sii]/Ha values of the order 0.20 ±0.02 (cf. Table 0) are 
also not predicted by this geometry. 

Additional information on ionization sources and re- 
sulting line ratios can only be made by using the simplified 
Ma86 model. Stars of spectral type 05 (T* = 4.5 • lO'' K) 
are taken to be responsible for the ionization of the ob- 
served elements. Taking error bars into account halo val- 
ues for [Nil] /Ha are well reproduced, as long as they are 
< 0.62. Larger ratios lead to a failure of the model. Rela- 
tive line strenghts of [S ii]/Ha (si and s2) are well fitted by 
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Fig. 14. NGC4634: Line ratios of [S ii] A6717/ Ha and 
[Nil] A6583/HQ; for slit s2. 5" correspond to 465 pc. The 
lower panel shows again [Sii]/[Nii] vs. EMhc along with 
the representative errors. 



Ma86. It produces values of 0.62 which is consistent with 
the measured data especially for slit position si. Although 
precise statements on excitation mechanisms of the extra- 
planar DIG cannot be made with the present data, the 
failure of the photoionization models implies more than 
one ionizing/heating source. 

4-7. Discussion 

Recent studies of DIG in the Milky Way and several edge- 
on galaxies confirm a nearly constan t tren d for [Sii]/Ha 
vs. [Nil] /Ha (e.g., H affner et al. 1999| ) and also for 
[Sll]/[Nii] vs. z (Rand 1998 ). These observations cannot 



be reproduced by photoionization models which depend 
on the ionization parameter U, because as [Nil] /Ha and 
[Sil]/Ha increase towards the halo, due to a diluted radi- 
ation field (smaller U), [S ii] / [N i i] inc reases, too. To ex- 
plain this finding Reynolds et al. ( |l999 ) recently proposed 
an additional heating source that is proportional to 
and increases the electron temperature towards the outer 
halo. 

Ideas for possible heating processes reach from tur- 
bulent dissipation to magnetic reconnection. In the low 
density environment these processes would be more effi- 
cient than photoionization, thus providing the necessary 
heating rate to account for the shapes of the above men- 
tioned line ratios and possibly also for the rise of [O iii]/H/3 
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Fig. 15. [Sii] A67I7/Ha vs. [Nil] A6583/Ha for all galax- 
ies with nitrogen and sulfur emission. The mean error for 
the halo of all galaxies is plotted in the upper right cor- 
ner . For comparison, data for the M ilky Way (Haffner et 
al. |1999[ ) and NGC891 (Rand |l998| ) would be located be- 
tween the dashed lines. 



~r" 



~r" 



" = NGC 1963 

■o = NCC3044 

■ = NGC440S 

♦ = NGC4634 



J I I I I I I I I I I I I I I I I I I I I I I I I I L 

-SO -lO O 10 20 



z (arcsec) 



Fig. 16. Alternative plot of [S ii] A6717/[N ii] A6583 vs. 
distance z perpendicular to the plane. Now, different gra- 
dients are visible in both halo hemispheres. 



with z (Reynolds et al. |I99E ). We present in Fig. 15 plots 
of [Six] /Ha vs. [Nii]/Ha supporting a linear dependence. 
Since line ratios are given with respect to Ha intensities, 
this relation holds with respect to densities. If [Sii]/[Nll] 
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line ratios are instead plotted vs. the geometrical unit z 
(Fig. 16) the scatter is much increased. In this presen- 
tation significant non-linear changes in [Sll]/[Nll] vs. z 
can be noticed. In some cases (NGC 3044, NGC 1963 and 
NGC 4634) the halo hemispheres may have different slopes 
which could be explained by, e.g., different metallicities 



(Haffner et al. 1999). Plots of line ratios vs. emission mea- 
sure EMffQ as a first estimate of a density dependence 
(lower panels of Figs. 4,9,12, and 14) also support a non- 
linear trend. 

As can be seen from Fig. 16 the ratio of [Sii]/[Nii] 
for NGC 1963 increases rapidly from 1.0 in the disk at 
z = 0" to 1.6 at z = —18" towards the halo. Values of 
1.12-1.34 are predicted by DM94 assuming the same pa- 
rameter setting as mentioned in Sect. 4.1 reflecting the 
observed data for the disk and the northern halo region 
well. For ratios less than 1 the chosen parameter set which 
previously reproduced [Nil] /Ha and [Sli]/Ha now leads 
to a model failure (southern halo, cf. Fig. 16). Alternative 
settings are also unable to fit observations. A likely phys- 
ical reason could be a small scale fluctuation (increase) of 
the electron temperature Te due to a decrease in metallic- 
ity or density of the extraplanar DIG. Although the model 
cannot account for local density or metallicity variations 
the general trend of [N ii] /Ha or [S ii] / [N ii] is qualitatively 
well reproduced. 

Compared to NGC 1963 the shape of [Sii]/[Nii] for 
NGC 3044 looks similar but reveals a lower starting value 
(0.67) and a non-linear gradient for \z\ < —10" (com- 
pare to Fig. 15 where the gradient seems to be linear). 
With respect to the mean errors, the observed ratios are 
again best represented by the DM94 model. The scatter in 
[Sii]/[Nii] (Fig. 16) at -20" < z < -10" suggests that the 
electron temperature and most likely also the metallicity 
is a function of z, varying significantly on large scales (« 
1 kpc). Within our sample NGC 1963 and NGC 3044 are 
showing the steepest gradients and largests mean values 
of [Sii]/[Nii], hence deviating from Rand's observation in 
NGC 891. 

The ratio of [Sii]/[Nii] for NGC4402 (Fig. 15 or 16) 
reveals only a relatively flat gradient towards the halo. 
This result differs extremely from the one for NGC 1963 or 
NGC 3044, and is similar to the ratio found for NGC 891. 
It also indicates that none of our selected models can pre- 
dict correct values, especially not for the most extreme 
points in this plot. It could be argued that diluted and 
very soft radiation is responsible for the observed values. 
However, this is not probable since [Nil] should be in- 
creased, too. 

Haffner et al. ( |1999D have shown for the Milky Way 
that different slopes in [Sii]/Ha vs. [Nil] /Ha could be an 
indicator for different metallicities. As a result the indi- 
vidual gradients visible in Fig. 16 could indicate differ- 
ent metallicities for each halo hemisphere. If this is true 
NGC 4402 would have the same metallicity as the Milky 
Way or its "twin" NGC 891. 



The data for NGC 4634 in Fig. 16 are in agreement 
with Ma86 and demonstrate a constant run at 0.6 for the 
disk, followed by a relatively steep gradient. At z = ±10" 
the ratio appears to remain constant at 1.0. In the disk 
[Sii]/[Nii] is comparable to that of NGC 891 but with 
increasing \z\ this similarity fades. 

In summary we note that strong nitrogen lines (> 0.62 
with respect to hydrogen) are generally not well repro- 
duced by pure photoionization models. The general trend 
of [Sii]/[Nli] is reproduced qualitatively but not quan- 
titatively (models cannot account for localized changes). 
They even fail in fitting oxygen and helium line ratios for 
the halo. The diagnostic diagram using [OIII] and [01] for 
NGC 1963 indicates that shocks may contribute as an ion- 
izing and heating source in the halo. This is similar to the 
finding of Martin ( |1997D for diffuse ionized gas in the out- 
fiows of dwarf galaxies. However, for the objects studied 
here, the discussed line ratios do not only change with ge- 
ometrical distance from the midplane of the disk, i.e. with 
respect to the location of the suspected ionizing sources. 
We rather find sudden localized changes, most evidently in 
NGC 3044 and NGC 4402, indicating that ionization con- 
ditions and/or heating rates may change on length scales 
of a few hundred parsecs. These non-linear changes can 
also be seen in the correlation between [SII]/[NII] line ra- 
tios and emission measure (Figs. 4,9,12, or 14). Since they 
also correlate with the local density of the diffuse medium, 
these variations can only be explained by small scale den- 
sity fluctuations. 

4-8. Kinematics 

In order to probe our results and to establish a relationship 
between excitation mechanisms and gas kinematics we 
have plotted in Fig. 17 an additional diagnostic diagram 
regarding line width (FWHM) vs. line ratio following, e.g., 
Lehnert & Heckman ( |1996| ). We have chosen [N ii] A6583 as 
kinematical tracer since its higher mass minimizes ther- 
mal broadening. Only three target galaxies with nitrogen 
emission of > 2 x 10~^^ erg s~^ cm~^ have been used 
for this plot. We point out that we did not correct for 
effects of kinematical line broadening. Thus, the data in 
Fig. 17 do not show the pure thermal line broadening. 
In order to establish nevertheless a relationship between 
ionization and gas dynamics we minimized the effect of 
kinematical line broadening in using only slit positions lo- 
cated near the center of each galaxy (si in NGC 1963, si 
in NGC 3044 and s2 in NGC 4634, cf. Fig. 0). A correction 
for the instrumental profile has been applied. 

The spectral resolution has been determined to be 4.6A 
by measuring the FWHM at different positions of HeAr 
calibration spectra. This translates into 210 kms~^ for the 
Ha emission line, corresponding to 2.3 on a logarithmic 
scale. Because most data in Fig. 17 are located between 
2.2 and 2.4 one has to be aware that only values larger 
than 2.3 (within the errors) have a significant physical 
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meaning. We therefore checked for consistency in measur- 
ing FWHM of the [Oi]A6300 nightsky line in NGC 1963 
and NGC 3044. Since this prominent emission feature is 
nearest to the [N ii] A6583 spectral line one thus obtains 
more precise values of the kinematical resolution. In the 
present case the parameter has to be modified slightly 
to 4.2 A or 200kms~^, leading again to the same value 
as mentioned above. Fortunately we are interested in line 
broadening towards the halo and hence the main body of 
the data can be used for analysis. 

For all galaxies plotted in Fig. 17 a clear correlation 
between gas kinematics and its ionization is visible. The 
broadening of lines at high galactic latitudes, even if less 
significant, supports a prominent large scale motion or in- 
creased turbulence of the ionized gas particles. Since this 
motion is tightly correlated to the ionization of the DIG 
it seems that in the halo of these galaxies shocks, super- 
winds or turbulent mixing layers as additional ionization 
or heating sources occur. 

This can be seen very nicely in NGC 1963 where disk 
and halo data cover different regions. Additionally, the 
hybrid character of NGC 3044 becomes obvious, again. 
Line broadening suggests only very moderate shock con- 
tributions. Although for NGC 4634 the data allows no de- 
tailed diagnostics concerning possible ionization mecha- 
nisms Fig. 17 suggests that shocks (besides photoioniza- 
tion) could represent one possible additional excitation 
mechanism of the extraplanar DIG. 



5. Summary and conclusions 

Line ratios for the halos of NGC 1963, NGC 3044 and 
NGC 4634 can be explained as a combination of pho- 
toionization by O stars and shock ionization. Since the 
measured FWHM increases with height above the galac- 
tic plane the results obtained from Fig. 17 support our 
findings that shocks can contribute to the ionization of 
the extraplanar DIG. The photoionization models from 
Mathis and Domgorgen & Mathis are not always capable 
to reproduce all fine ratios correctly (high [N ii] A6583/Ha, 
[Oi] A6300/Ha, and [Hei] A5876/Ha). These models were 
developed to describe the measured line ratios for the 
Milky Way (Reynolds 19854 b). However the local envi- 
ronment near the sun is not necessarily representative for 
the more active regions in our target galaxies. Thus new 
models are required which should rest upon recently pub- 
lished extragalactic data and consider additional ioniza- 
tion/heating mechanisms such as shocks. 

It is important to mention that the ratios for all galax- 
ies in Fig. 16 are not constant over the whole range of 
z and can vary significantly on small scales (w 250 pc). 
These results deviate from observations made e.g., by 
GoUa et al. ( |l996[ ) Ran d ( |199^ ), Otte& Dettmar ( [1999D , 
or Haffner et al. (1999). Photoionization models can ex- 
plain the general shape of [Sii]/[Nii] but cannot account 
for individual changes (e.g., NGC 4402). In addition non- 
linear changes, as visible in Fig. 16, correlate with density 
and are explainable by small scale density fluctuations. 

If variations in S"*" reflect variations in electron tem- 
perature, leading to = T^^z), these fluctuations could 
be due to changing metal abundances or varying densities 
within the DIG. If different slopes in Fig. 15 or 16 are 
indicators of different metallicities (Haffner et al. 1999), 
NGC 4402 would have the same metallicity as the Milky 
Way or NGC 891. In accordance with Lehnert & Heckman 
( 1996| ) a correlation between line ratios (ionization) and 
line widths (gas kinematics) could be established. Our re- 
sults of a decrease of [Oili]/H/3 with increasing \z\ and 
a significant scatter in [Sii]/[Nii] are different compared 
to recent measurements in NGC 891. All these aspects in- 
dicate that the ionization of the extraplanar DIG varies 
from galaxy to galaxy and within a galaxy. 
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